A numerical finite difference method is developed to simulate the viscous flow over blunt/sharp bodies at incidence. Herein, a two-layer model is suggested. The inner region consists of the three-dimensional boundary layer and boundary region. The governing system applies in boundary regions and for problems with cross flow reversal. The equations are integrated by a predictor-corrector scheme. For the turbulent boundary layer analysis, both a mixing length model and a two-equation kinetic energy-dissipation system is considered for Reynolds stress closure. Numerical solutions are compared with experimental data and indicate that the present formulation can give an accurate prediction of aerodynamic loads, skin friction and heat transfer rates on sphere-cone-cylinder-flare shape bodies at angle of attack. The calculations are suitable for (i) supersonic or hypersonic freestreams, (ii) large Reynolds number and (iii) flows without streamwise flow separation; however, secondary flow reversal is allowed.